Analytic determination of dynamical and mosaic length scales in a Kac glass model 
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We consider a disordered spin model with multi-spin interactions undergoing a glass transition. 
We introduce a dynamic and a static length scales and compute them in the Kac limit (long-but- 
finite range interactions). They diverge at the dynamic and static phase transition with exponents 
(respectively) —1/4 and —1. The two length scales are approximately equal well above the mode 
coupling transition. Their discrepancy increases rapidly as this transition is approached. We argue 
that this signals a crossover from mode coupling to activated dynamics. 

PACS numbers: 05.20.-y (Classical statistical mechanics), 64.70.Pf (Glass transitions), 75.10Nr (Spin-glass 
and other random models) 



The relaxation mechanisms of supercooled liquids near 
the glass transition are today poorly understood. For 
moderate supercooling, mode coupling theory (MCT) 
captures several important dynamical features, such as 
time scale separation and the relation between a and (3 
relaxation Unfortunately MCT fails at low temper- 
ature since it predicts a spurious power law divergence 
of the relaxation time at an crgodicity breaking temper- 
ature Tfj. 

A partial elucidation of ergodicity restoration below 
Td originated from the remark that MCT is exact for 
some mean- field disordered spin models In this con- 
text, the relaxation time divergence at Td finds its root 
in the proliferation of metastable states with diverging 
lifetime. The existence of such states provides the ba- 
sic ergodicity restoration force in finite dimension. Phe- 
nomenological droplet arguments based on this mean 
field picture 0, 0] predict in fact the existence of about 
explE^''} different states at length scale £, while a given 
boundary condition can select one of them by raising the 
free energy of the others by at most = al^ with 
< d — 1. Balancing these two forces yields a "mosaic 
length" ~ (/^f/S)^^, above which the system is ther- 
modynamically a liquid. 

The above picture suggests that the typical size of 
ergodicity restoring rearrangements is determined by 
thermodynamics. According to an alternative intuition 
(largely inspired by kinetically constrained models @), as 
the temperature is lowered, it becomes impossible to re- 
lax local degrees of freedom through purely local moves. 
The relevant length scale is essentially dynamical and 
corresponds to the size of cooperative rearrangements. 

In both scenarios, the typical size of rearranging re- 
gions can be characterized by a dynamical (four points) 
susceptibility and a corresponding dynamical length 
@, 0] ■ Recent experiments on colloidal systems and su- 
percooled liquids have indeed revealed a marked increase 
of ^d close to the glass transition m. Within MCT, ^d 
has been shown to diverge at Td [a, ^ with no static 



counterpart. 

A relation between mosaic length and the dynamical 
length ^d was first noticed by Jack and Garrahan in 
a plaquette glass model. However a detailed theoretical 
description of this relation is lacking. 

The aim of this letter is to present the first analytic 
calculation of the mosaic length scale, and its compar- 
ison with the dynamical length computed in a unified 
framework. To this end we study a disordered model 
with Kac interactions [l^, and define length scales 
through properly chosen point-to-set correlation func- 
tions. We are thus able to answer some fundamental 
open questions, such as their relative order of magnitude 
in different temperature regimes. 

It is well known that conventional static correlation 
functions do not show any signature of a large length 
scale in supercooled liquids. In order to circumvent this 
problem, we define length scales through "point-to-set" 
correlations 0, We fix a reference configuration 

= {si"^}a;gz<i drawn from the equilibrium Boltzmann 
measure, and considers a second configuration s that is 
constrained to be close to s^"^ outside a sphere B^g {£) of 
radius £ centered at a particular point xq . In introducing 
the static correlation length is, s is distributed accord- 
ing to the Boltzmann law inside Bxq (i) with a boundary 
condition determined by sj-"\ We then define is to be 
the smallest value of i such that the correlation between 
s^"-* and s at the center of B^g (i) (as measured through 
a correlation function, e.g. (sio''s2;o)conn) decays below a 
preassigned value e. 

Consider now a particular choice of local stochastic 
dynamics verifying detailed balance. Let r be the auto- 
correlation time for a local degree of freedom. As shown 
in Ref. in a finite range system, one has necessar- 
ily Ci4 < T < exp(C2^f), with Ci and C2 two model 
dependent constants. The lower bound corresponds to 
a MCT-likc dynamics, a good description for T 3> Td 
(up to a non-trivial scaling exponent, i.e. r ~ £| with 
z > 1). The upper bound corresponds instead to the ac- 



2 



150 



120 



90 



60 



30 











- 






'"-^^^ - 

J , , : 




■ 




















— ■ — 1 











0.75 



0.78 



0.81 



0.84 



FIG. 1; Dynamic (empty squares) and static (filled squares) 
rescaled lengths Id and Is, for a Kac model with two and four 
spin interactions. The vertical lines correspond to the mode 
coupling and thermodynamic glass transitions Td « 0.813526 
and Ts ~ 0.766287. In the inset, the rescaled lengths as func- 
tions of (T - Td) (for ?d) and (T - T^) (for 4). The dotted 
lines have slope (respectively) —1/4 and —1. 



tivated behavior that should hold below Td (we expect 
T ~ exp(^^) with ^ < d — 1). In a model with large- 
but- finite interaction range (see below), one finds 
£s{T,j) « 7~^€s(r), and the two types of dynamics can 
be distinguished sharply (exponentially in I/7). 

The two regimes can be bridged defining a new dy- 
namical length id which separates scales £ ^ on 
which non-activated relaxation is possible, from scales 
£ <C ^d on which activation is necessary. In order to pre- 
cise this notion, consider a system initialized by setting 
s(t = 0) = s*-"-* and constrained to remain close to s^"^ 
outside {£) at all subsequent times, and let t{£) be the 
shortest time such that {sxo iO)sxo (O)conn decays below e. 
By the above definition t{1) =00 for £ < ig. We define 
id by the property that t{£) is (exponentially) divergent 
as 7 ^ for £ < £d, and stays bounded for £ > i^. 

In the following we shall compute ig and £d in a d- 
dimensional spherical p-spin disordered model with Kac 
interactions. Its elementary degrees of freedom are soft 
spins Si G M, associated to the vertices of a cubic lattice of 
size 2L (i.e. i e {—L + 1, . . . , L}'^), with periodic bound- 
ary conditions. Given the interaction range 7"^ > and 
a non-negative rapidly decreasing function ip : R_|_ — > , 
normalized by /'(/'(Nl) d'^x = 1, we define the local 
overlap of two configurations sj-^^ and sj-^"^ as Qi2{i) = 



(1)J2) 



The local spherical constraint 



is Qn(«) = 7'^Ej V'(7N - j\) = 1- Hamilto- 

nian H{s) is a quenched Gaussian random variable of zero 
mean and covariance E[i/(s(i))iJ(s(2))] = Y^j{Qi2{i)) 
where f{x) is a polynomial with positive coefficients. In 
the following we shall use as running example the poly- 
nomial f{x) = x"^ +x'^ (i.e. a model with two and four 
spin interactions) in d = 2 dimensions (but the critical 



behavior is independent oi d > 2). While the quartic 
term in the polynomial assures the wanted phenomenol- 
ogy at the mean-field level, the relatively small quadratic 
term is introduced to have a regular gradient expansion of 
the free-energy functional (see below and |l11|]). The Kac 
limit is defined as T ^ 7^^ » 1 (in other words, the limit 
7 ^ is taken after L 00). It is convenient to define 
the rescaled lengths -^g/d = 7 -^s/d that admit a finite limit 
as 7 ^ 0. In Fig. [T] we plot £d and 4 as functions of tem- 
perature for our running example, in the Kac limit. They 
diverge at two distinct temperatures Td and Tg. The dy- 
namical length diverges as £d ~ (T — Td)"^/*. This is the 
same behavior as (independently) found within MOT for 
the length ^d @, we therefore identify id with ^d- 
The static length behaves as 4 ~ 1/S(T) ~ (T - T^)"! 
corresponding to 9 = d ~ 1. 

We shall now explain how the lengths id and is have 
been computed. Consider a reference configuration s^°'\ 
a sphere Bo{£) with radius £, centered at x = 0, and let 
q < 1. We introduce the constrained Boltzmann measure 
{■)" by setting 

= ^ j 0{s) e-P"^^'> n '5(Q.,.(=)W-g), (1) 

where denotes integration over configurations satisfy- 
ing the local spherical constraint. For technical reasons 
it is preferable to take q < 1: s is only required to be 
close to s*^"^ outside Bo(f). The results are qualitatively 
independent of q if this is large enough. Next we define 
the correlation function 



G(£)^E{E,<„, [(Q,^,(c,(0))?]} 



(2) 



where Eg(Q) and E denote (respectively) averages over the 
reference configuration and the quenched disorder. The 
static length scale £s is the smallest i such that G{i) < e. 

The averages Eg(Q) and E can be taken by introducing 
replicas along the lines of Uj. Integrals over the spin 
variables are then traded for and nr x nr matrix order 
parameter qab{i)- Next we rescale the position to define 
X — ji ^ [— i, L]'^, L = 7L, and write (with an abuse of 
notation) qat = qab{x), to get 



Gii) 



lim 



9i,n+i(0) e 



d[qab 



(3) 



The dependency upon 7 is now completely explicit and 
the functional integral can be performed using the sad- 
dle point method. Inspired by the solution of the mean 
field model 15|, we look for a one-step replica sym- 
metry breaking (IRSB) saddle point ql^^^{x). This 
is characterized by three scalar functions p{x), qi{x) 
and qo{x) and by a single Parisi IRSB parameter m. 
While p{x) is the local overlap between the reference 
configuration s*-"^ and the constrained one s, qi{x) and 
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FIG. 2: Saddle points of the IRSB action, cf. Eq. © for 
T — 0.81365 and several values of the ball radius £. Here 
we plot the outcome of the iterative solution initiated with 
the p{x) = q identically. For £ < = 59, this is the high- 
overlap solution. For > no high-overlap solution exists; 
the remaining dependence upon £ is essentially a shift. 

Qoix) are the overlaps of two constrained configurations 
when they belong (respectively) to the same or to dif- 
ferent metastable states. Using this ansatz, one obtains 
S[Qab] = nj C{x) d'^x + 0{n^), where 



C = 



{/(I) + 2f{^ *p)-{l- m)fi^P * gi) 



P - 90 



2 1 — (1 — 'm)qi — mqo 



(4) 



m — 1 1 

log(l - qi)- — log 1 - (1 - mjqi - mqo\ , 

2m Zm 

with the various fields being evaluated at position x. Here 
we used the shorthand ip * q{x) — J ip{\x — y\) q{y) d'^y. 
The constraint enforcing s to be close to s*^"^ outside 
Bo(^) is fulfilled by setting p{x) —qioix^ Bo(^) (notice 
the rescaling of the radius: i — •ji). 

Using Eq. ([3]), we obtain liuiy^Q G{£/j) = p{0), where 
p{x) is evaluated at the minimum of the action 5*0 = 
JC{x)d'^x. This imphes that, as anticipated, £s{T) 
"f£s{T, 7) has a finite 7 — > limit. We minimize the action 
looking spherically symmetric saddle points and iterating 
numerically the Euler-Lagrange equations for ^ until a 
stationary point is reached. We then repeat for several 
initial values of p{x), qo{x) and qi{x) conditions, and 
compare the corresponding stationary points. 

Since this procedure is relatively heavy from a com- 
putational point of view, we simplify the Lagrangian 
(jU in two ways. First, we expand the terms of the 
form * q)ix) in gradients of q{x), and truncate to 
second order, thus obtaining f{q) + cf'{q)Aq{x), where 
^ ~ Id J d'^z (in our running example c = 1). 
Second, we discretize the model on a lattice of spacing 
h in the radial coordinate. Neither of these simplifica- 
tions is expected to modify the qualitative behavior of 



FIG. 3: Plot of the overlap p{0) as a function 
of £ in the approximated model with for T = 
1.0101,0.9065,0.8319,0.8222,0.81510,0.8142 > Td and T = 
0.7990 < Td. Inset: comparison between the approximate 
(analytical) and complete (numerical) minimization of the ac- 
tion for T = 0.8222. 

the model, and we checked our predictions against alter- 
native discretizations. 

In Fig. [2] we show some typical results for T > T^. 
For small £, a unique solution is found. The overlap p{x) 
between the probe replica and the reference one is every- 
where large: close to q on the border, it saturates well 
inside the ball to a smaller value gd roughly independent 
of the ball radius. The two overlaps qoi^x), qi{x) are equal 
and close to p{x). As £ crosses some small value £0, a new 
solution appears, with p{x) rapidly decaying to in the 
interior of the ball, and the influence of the boundaries 
remaining confined to a small distance p. This solution 
describes a constrained replica s well decorrelated from 
the reference one s*-"' , but has initially higher value of 
the free energy and is therefore thcrmodynamically ir- 
relevant. The static length £s is the smallest value of £, 
such that the free energy of the low-overlap solution be- 
comes smaller than the free energy of the high-overlap 
one. Above this length, the correlation between s and 
_ gioi) in the center of the ball is small. However, the high- 
overlap solution persists as a metastable state. If the 
system is let evolve from the initial condition s^") (and is 
constrained outside the ball), it will take an exponential 
(in 1/7) time to equilibrate. The high-overlap solution 
disappears at a larger radius £d- For £ > £d, the system 
is no more trapped in a metastable state selected by the 
reference configuration. 

The evolution of different solutions as a function of 
£ can be followed on Fig. [3l where we plot the overlap 
p(0) between s^"^ and s at the center of the ball. For 
T ^ Td, p{0) is a single valued function of £. Closer 
to Td, multiple branches develop, with the coexistence 
region diverging as T | Td. Beyond the high- and low- 
overlap solutions described above, an intermediate un- 
stable branch is also present. In order to recover the full 
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FIG. 4; Typical shape of the configurational entropy 
for 4 < ? < ?d (here T = 0.814 > Td, 4 ~ 16, id 41, 
and £ — 25). Inset: configurational entropy E^(/*)/£'* of the 
dominating metastable states as a function of i. The curve 
crosses the axis at l'-^'^^ and ends at l^- 

curve (including the unstable branch) we resorted to a 
simple approximation: we minimized the action under 
the constraint p{x) — qo{x) — qi{x). In this setting, and 
under the "thin wall approximation" yjj] , the problem is 
equivalent to a one-dimcnsional mechanical system, and 
can be solved by quadratures. As shown in the inset of 
Fig. [3l the approximate and complete actions give results 
in good agreement. In particular, both cases yield the 
critical behaviors 4 - (T-Ts)"! and £d (T-Ts)-^/^. 
Below Td all three branches extend to infinity, indicating 
that activation is necessary on all length scales. 

The complete IRSB lagrangian ^ allows for a more 
precise description of the system. Consider T close to Td 
and id [£o, id] where io is the minimal length for multiple 
solutions. As anticipated, on the low overlap solution, we 
find gi > go- The interpretation is, as usual, that this 
solution does not describe a liquid phase, but rather a 
glassy one, with many pure states. The internal overlap 
of each such state is qi and the overlap between distinct 
states go- The number of metastable states is exponential 
in £ — i/j for £ > i^^^^ (the action is maximized by 
m — 1) and subexponential for £ < fii^SB (^j^, < i) j^gg 
slightly modify the free energy of the low overlap solution 
for i < ^ir'SB -(-j^^g changing the value of is- 

Let us reconsider the thought experiment of letting the 
system evolve starting from the reference configuration 
gia) ^ ^j^g system escapes the metastable state selected by 
s^"-* in a time exponential in I/7. After that, however, it 
does not equilibrate freely within the rest of phase space, 
but rather moves through activation from one state to 
the other, spending an exponential time in each one. 

By taking the Legendre transform of the m-dependent 
free energy, we can compute the configurational entropy 
E^(/) (i.e. the log number of metastable states as a func- 
tion of their internal free-energy /) [? ]. In Fig. [?] we 



show an example of its behavior for £s < i < id- An 
exponential number of metastable states is present in an 
^dependent free energy interval [fi , (most of them 
having free energy /2). The dominating ones have free 
energy determined by maximizing S^/) — /?/• The 
free energy T+ of the high overlap state is lower than for 
most of other metastable states < Ji (because of the 
better matching with the boundaries) but higher than the 
overall free energy T_ = /» — TS^/*). If i is decreased, 
the gap between F- and T^ decreases. At ^u^SB ^^^^ 
figurational entropy E^(/,) vanishes and F- = /* ~ fi- 
Eventually T_ and F+ cross at ?= 4 < ?iRSB_ ^^le 
other hand i is increased, the configurational entropy in- 
creases as well, until the IRSB solution ceases to exist. 
This happens at the dynamical length i^-, coherently with 
the fact that beyond this scale relaxation does not require 
activation. 

For a small non- vanishing 7, the above picture should 
remain roughly correct. However, activated and MCT 
time scales must now be compared. The first one corre- 
sponds a typical relaxation time Tact ^ exp{T£f }, while 
the second gives tmct ~ ^d- Therefore, MCT behavior 
is faster and dominates until very close to Td. The di- 
vergence is then cut-off when id{T) « exp{T4(Td)'^/2;}, 
and dynamics becomes activated at lower temperatures. 
During this cross-over, the physical dynamical length 
characterizing the size of cooperatively rearranging re- 
gions crosses over from ~ ^d, to ~ £s- 

Summarizing we implemented the mean-field theory 
of mosaic state in a microscopic model. This allowed to 
establish relations between mosaic and MCT dynamical 
lengths which had not been predicted from phenomeno- 
logical considerations. The dynamical length scale is de- 
fined by the property that dynamics is dominated by ac- 
tivation at shorter scales. It can further be identified 
with the one appearing in the MCT divergence of four 
point susceptibilities. Below mosaic length the system is 
instead rigid for thermodynamic reasons. The two hap- 
pen to be close above the MCT critical temperature and 
widely separated between this and the glass transition, 
corresponding to the crossover between relaxational and 
activated dynamics. 
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